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Modification of the right-hand-side of canonical commutation relations (CCR) naturally occurs if 
one considers a harmonic oscillator with indefinite frequency. Quantization of electromagnetic field 
by means of such a non-CCR algebra naturally removes the infinite energy of vacuum but still results 
in a theory which is very similar to quantum electrodynamics. An analysis of perturbation theory 
shows that the non-canonical theory has an automatically built-in cut-off but requires charge/mass 
renormalization already at the nonrelativistic level. A simple rule allowing to compare perturbative 
predictions of canonical and non-canonical theories is given. The notion of a unique vacuum state is 
replaced by a set of different vacua. Multi-photon states are defined in the standard way but depend 
on the choice of vacuum. Making a simplified choice of the vacuum state we estimate corrections to 
atomic lifetimes, probabilities of multiphoton spontaneous and stimulated emission, and the Planck 
law. The results are practically identical to the standard ones. Two different candidates for a 
free-field Hamiltonian are compared. 



I. INTRODUCTION 



The standard quantization of a harmonic oscillator is based on quantization of p and q but w is a parameter. To 
have, say, two different frequencies one has to consider two independent oscillators. On the other hand, it is evident 
that there can exist oscillators which are in a quantum superposition of different frequencies. The example is an 
oscillator wave packet associated with distribution of ccnter-of-mass momenta. It is known that the superposition of 
momenta gets translated into a superposition od Doppler shifts and therefore also of frequencies. We stress here the 
word "quantum" since the superpositions we have in mind are not those we know from classical oscillations. 

This trivial observation raises the question of the role of superpositions of frequencies for a description of a single 
harmonic oscillator. The motivation behind the problem is associated with the question of field quantization: Is it 
possible that a quantum field consists of oscillators whose frequencies are indefinite? If so, maybe to quantize the 
field it is sufficient to use only one oscillator which exists in a quantum superposition of all the possible frequencies 
allowed by the boundary conditions of a given problem? 

The idea is very simple. It is known that a "one-particle" state vector can be regarded as a representation of an 
ensemble of particles in a given pure state. On the other hand, the classical electromagnetic field can be regarded 
as an ensemble of oscillators. The standard idea of quantization, going back to 1925 is to treat the field as an 
ensemble of quantum oscillators. But the ensemble itself is, in a sense, a classical one since for each frequency we need 
a separate oscillator. This is analogous to a classical ensemble of particles forming a classical wave on a lake surface. 
For each point on the surface we need a separate particle because a classical particle can ocupy only a single point in 
space. A quantum wave is of course different and we are all accustomed to the idea of a single-particle wave. In this 
case the properties of the entire ensemble are somehow encoded in properties of a single element of the ensemble. 

For some reasons, probably partly historical, it seems that the idea of a single-particle state vector representation 
of the ensemble of oscillators has never been considered. The historical reason may be the fact that the very concept 
of field quantization occured already in 1925. At that stage quantum mechanics existed still in a matrix form and 
the Schrodinger paper "Quantisierung als Eigenwertproblem" where the Schrodinger equation occured for the 
first time and the role of eigenvalues was explained, was not yet published. Actually, as explained by Jammer ||, 
Heisenberg quantized the harmonic oscillator without having heard of matrices so no wonder he could not treat the 
parameter u> entering E = p 2 /2m + muj 2 q 2 /2 as an eigenvalue of some operator. 

Heisenberg's quantization leads to the well known algebra of canonical commutation relations (CCR) [a w , ajj = 1, 
where 1 is the identity operator and lo a classical parameter. As it turns out the replacemwnt of u> by an operator Co 
leads to non-canonical commutation relations (non-CCR) [a w , a/J\ = l w , where u> is an eigenvalue of Co and ^ w l w = 1. 
l u , similarly to 1 commutes with all creation and annihilation operators and the remaining commutators of non-CCR 
are the same as those of CCR. This subtle difference of the right-hand sides of CCR and non-CCR immediately 
explains why a non-CCR-quantized electromagnetic field will have vacuum with finite energy. Since the non-CCR 
and CCR algebras are so similar to each other it is not surprizing that the resulting theories are also very similar. 
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The main consequence of the non-CCR modification of CCR is different normalization of n-photon states. For 
example, 

£(010*410) = 1 (1) 

and therefore (Ola^aJjO) < 1. This should be contrasted with the CCR result (Ola^aJjO) = 1, and the resulting 
divergence 

£<Qk,4 10) = oo (2) 

UJ 

which is the source of infinite vacuum energy. 

To end these introductory remarks one should mention that several approaches towards an alternative description of 
the electromagnetic field at a fundamental level were already proposed (e.g. Janes' j| neoclassical theory, stochastic 
electrodynamics || ) . But the main idea of all such alternatives was to treat the field in classical terms and to associate 
the observed discreteness of emission/absorbtion phenomena with the quantum nature of atoms and not with the field 
itself. 

The approach we will discuss in this paper does not belong to this tradition, is much more radical and, so to say, 
goes in the opposite direction. We will not try to make the field more classical. What we will try to do is to make it 
even more quantum by replacing classical parameters with eigenvalues. 



II. HARMONIC OSCILLATOR IN SUPERPOSITION OF FREQUENCIES 

We know that frequency is typically associated with an eigenvalue of some Hamiltonian or, which is basically the 
same, with boundary conditions. A natural way of incorporating different frequencies into a single harmonic oscillator 
is by means of the frequency operator 

q = u k \uj k ,j k )(uj k j k \ (3) 

u>k,3k 

where all uj k > 0. For simplicity we have limited the discussion to the discrete spectrum but it is useful to include 
from the outset the possibility of degeneracies, represented here by the additional discrete quantum numbers j k . The 
corresponding Hamiltonian is defined by 

H = htt®^(a^a + fto') (4) 



where a — J2^Lo + l\n)(n + 1|. The eigenstates of H are \uj k ,j k , n) and satisfy the required formula 

H\uJ k ,j k ,n) = Huj k (n+ ^j\u; k ,j k ,n) (5) 

justifying our choice of H . The standard case of the oscillator whose frequency is just u> coresponds either to = wl 
or to the subspace spanned by \ui kl j kl n) with fixed w k = u>. Introducing the operators 

a ^k,j k = \uk,jk}(ukjk \ <8>a (6) 

we find that 

H = \ J2 ^ (°-* Jk J* + A a L Jh ) • ( ? ) 

The algebra of the oscillator is "noncanonical" : 

[aukMiol^jJ = \uk,jk){uk,3k\ ® 1 = low* (8) 

K* Jk > a l t ,j, } = for H > 3k) ¥= (Wi ,3l) (9) 

[au )h ,j h ,a Ul j l ]=0 (10) 

[«L, J v^ ! ,iJ= Q ( n ) 
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The dynamics in the Schrodinger picture is given by 



ihdt\V) = H\V) = (a t a + ^l)|*>. (12) 

In the Heisenberg picture we obtain the important formula 

a^ k (t)=e iHt/R ^,j k e- iHt/R (13) 

= \u k , jk)(uk,jk\ ® e-^a = e- iUkt a Uktjk . (14) 

Taking a general state 

W)= ^2 ip{uk,jk,n)\uj k ,j k ,n) (15) 

we find that the average energy of the oscillator is 

(H) = (i/>\H\1>) = ^ |^ fe ,i fe ,n)| 2 na; fe (n+i). (16) 

The average clearly looks as an average energy of an ensemble of different and independent oscillators . The ground 
state of the ensemble, i.e. the one with tp(tuk,jk, n > 0) = has energy 



(H) = \ ]T |VH,j fc ,0)| 2 fiw fc (17) 

which is finite if 
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^2 ^(wfc, jfc,0)|w fc ,jfc) (18) 



belongs to the domain of tt. The result is not surprising but still quite remarkable if one thinks of the problem of 
field quantization. 

The very idea of quantizing the electromagnetic field, as put forward by Born, Heisenberg, Jordan [Q and Dirac || , 
is based on the observation that the mode decomposition of the electromagnetic energy is analogous to the energy of 
an ensemble of independent harmonic oscillators. In 1925, after the work of Heisenberg, it was clear what to do: One 
had to replace each classical oscillator by a quantum one. But since each oscillator had a definite frequency, to have 
an infinite number of different frequencies one needed an infinite number of oscillators. The price one payed for this 
assumption was the infinite energy of the electromagnetic vacuum. 

The infinity is regarded as an "easy" one since one can get rid of it by redefining the Hamiltonian and removing 
the infinite term. The result looks correct and many properties typical of a quantum harmonic oscillator are indeed 
observed in electromagnetic field. However, subtraction of infinite terms is in mathematics as forbidden as division by 
zero so to avoid evident absurdities one is forced to invent various ad hoc regularizations whose only justification is 
that otherwise the theory would not work. In larger perspective (say, in cosmology) it is not at all clear that an infinite 
(or arbitrarily cut off at the Planck scale) energy of the vacuum does not lead to contradictions with observational 
data 0. Finally, Dirac himself had never been fully satisfied by the theory he created. As Weinberg put it, Dirac's 
"demand for a completely finite theory is similar to a host of other aesthetic judgements that theoretical physicists 
always need to make" Q . 

The oscillator that can exist in superpositions of different frequencies is a natural candidate as a starting point for 
Dirac-type field quantization. Symbolically, if the Heisenberg quantization is p 2 + w 2 q 2 i— > p 2 + uj 2 q 2 , where w is a 
parameter, the new scheme is p 2 +ui 2 q 2 i— > p 2 + Cu 2 q 2 , where Cj is an operator. Its spectrum can be related to boundary 
conditions imposed on the fields. 

We do not need to remove the ground state energy since in the Hilbert space of physical states the correction is 
finite. The question we have to understand is whether one can obtain the well known quantum properties of the 
radiation field by this type of quantization. 
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III. "FIRST QUANTIZATION" — ONE-OSCILLATOR FIELD OPERATORS 

The new quantization will be performed in two steps. In this section we describe the first step, a kind of first 
quantization. In next sections we shall perform an analogue of second quantization which will lead to the final 
framework. 

The energy and momentum operators of the field are defined in analogy to H from the previous section 

H = k)(s, k\ \ (a) a + aa)^\ (19) 

s.ii 

= i ( a U a sM + asMala) ( 2 °) 

S.K 

P = ^M\s, k){s,k\ ^(^a + aa^ (21) 

S.K 

= \j2 M ( a U a sM + a s,Kal,n) ( 22 ) 

where s = ±1 corresponds to circular polarizations. Denote P = (H/c, P) and P ■ x = Ht — P ■ x. We employ the 
standard Dirac-type definitions for mode quantization in volume V 



l(t, x) = 2 (a^e-^e^e^ + <,e^'e s >-«-) (23) 

S.K ' 

_ e iP-x/hj£ e -iP-x/h ^24) 



E(t, x) = iJ2 ] /^- (a^e-^e^e^ - a^e*"^-'*^) (25) 

S.K 

= e iP^/^Ee- ip - x/n (26) 

(27) 

^ f ) = 1 E VW" K X {^-^^es,, - al^e-^el,) (28) 

S.K 

= e iP-x/*-~Be- ip - x/h , (29) 



where 



a s,fi — \s, k)(s, k\ ® a (30) 
a\ R = \s, k)(s, k\ <g> a f . (31) 

For later purposes we introduce the notation 

[a s ,z, a[-] = Ism = \s,k) (s, k\ <g> 1. (32) 

Now take a state (say, in the Heisenberg picture) 

I*) = E *s,z,n\s,K,n) (33) 

s,K,n 

= J2<S> s M\s,K)\a s ^) (34) 

where |a S) g) form a family of one-oscillator coherent states: 

a\a s ,a) = a s ^\a s ^) (35) 

The averages of the field operators are 
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<*|l(i,a!)|*) = E l^y^a^e— e s , R - <*e«'*e**) (37) 



(*| J(t,a!)|*) =*El^| 2 ^(^ K I « K x e s , K - - < s e— n K x e^) (38) 

These are just the classical fields. More precisely, the fields look like averages of monochromatic coherent states with 
probabilities |$ Si g| 2 . The energy-momentum operators satisfy also the standard relations 

H =\ J d 3 x^E(t,x)- E(t,x) + B(t,x)- B(t,xfj, (39) 

P= d 3 xE(t,x) x B(t,x). (40) 
Jv 



IV 

To end this section let us note that 



= E ^I<M 2 (KkI 2 + ^) (4i) 

S.K 

(*|P|*)=^M$^| 2 (K K H 2 + ^)- (42) 

The contribution from the vacuum fluctuations is nonzero but finite. One can phrase the latter property also as 

follows. The noncanonical algebra of creation- annihilation operators satisfies the resolution of identity 

wheras the canonical algebra would impliy 

EK^U^ 001 - ( 44 ) 



IV. "SECOND QUANTIZATION" 

The Hilbert space of states of the field we have constructed is spanned by vectors \s,K,n). Still there is no doubt 
that both in reality (and the standard formalism) there exist multiparticle entangled states such as those spanned by 
tensor products of the form 

|+,M,1)|-,k 2 ,1), (45) 

and the similar. It seems that there is no reason to limit our discussion to a single Hilbert space of a single oscillator. 
What we have done so far was a quantization of the electromagnetic field at the level of a "one-particle" Hilbert 
space. Similarly to quantization of other physical systems the next step is to consider many particles. What is not 
obvious (physically) is whether the oscillators should be considered as noninteracting. This physical freedom leads to 
two natural candidates for a free-field Hamiltonian. 

The noninteracting extension is essentially clear. Having the one-particle energy-momentum operators P a (i.e. 
generators of 4-translations in the 1-particle Hilbert space) we define in the standard way their extensions to the 
Fock-type space 

Va = Pa 

e(p a ®i + i(g)P a ) 

e(p a ®i®i + i®p a ®i + i®i®p Q ) 

© . . . . (46) 
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The x-dependence of fields is introduced similarly to the one-particle level 

F(t, x) = ei-P-x/Hf e -iV-x/h ( 4? ) 

but the field itself has yet to be defined. Assume 

T = ciF 

ffic 2 (F© 1 + 1 OF) 

0c 3 (F®1® 1 + 1®.F® 1 + 1® 1®P) 

© . . . (48) 

where Cfc are constants discussed below, and F is A, E, or _B. The multi-oscillator annihilation operator associated 
with such fields must be therefore of the form 

©c 2 (a s ,K (gil + 1® o a ,g) 

©c 3 (a S) g ® 1 ® 1 + 1 ® a s> g ® 1 + 1 ® 1 ® o a ,g) 

© . . . . (49) 

Having two 1-particle operators, say X and Y, one can easily establish a relation between the 1-particle commutator 
[X, Y] and the commutator of the extensions X, y-. 

[x,y]=c\[x,Y] 

®4([x, y]®i + i®[x, y}) 

©c| ([X, Y] ® 1 ® 1 + 1 ® [X, Y] ® i + 1 ® i ® [x, y]) 

© . . . . (50) 
The annihilation operators so defined satisfy therefore the algebra 

Ka,o^,]=0 for (s,k) # (s',k'). (5!) 

[a s ,g,a^] = l s ,g, (52) 

[a s ,B,a s ',g'] = (53) 

=0 (54) 

where the operator l Si g is defined by 

1 - - r 2 1 - 

©c|(l Si g ® 1 + 1 ® l s ,g) 

©<|(l s>S © 1 © 1 + 1 © l Sj g © 1 + 1 © 1 © l S) g) 

©..., (55) 

and l Si g is a single-oscillator operator (|32"|). 

An important property of the 1-oscillator description was the resolution of identity (43). The requirement that the 
same be valid at the multi oscillator level leads to c„ = 1/y/n. In such a case one finds that 

1* a + l.,a (56) 

but nevertheless 

$3l.,* = l, ( 5? ) 

that is l Sl g are the so-called positive operator valued (POV) measures ||. Below we shall give still another justification 
of this particular choice of c„. 
We can finally write 
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A(t, x) = J ^7 (a^e-^e^e^ + a^e^^e-***) 

S,K V 

— e iV-x/h^ e -iV-x/h 

.) 



B{t,x) = «E V x ~ "U^^e-^'^e^) 

_ e iV-x/hg e -iV-x/h 

These operators form a basis of the modified version of nonrclativistic quantum optics. 

Let us return for the moment to the case of a general c n . A straightforward calculation shows that 



H 



i J d 3 x(i(t,x)-£(t,x)+B(t,x)-B(t,x)J = ^^hug(al~a s ,z + a s ,zal^ 



© c 2 (h,{ a s,R, a\ R } + \{a s ,z, a\-} + a\ R ® a Si u + a s ^ ® a\ Sj 

© c 3 (H a «,«' a !,«} O 1 O 1 + 1 O H a «,«' a l,«} ®1 + 1®1® H a ^«' a l,«} 
+ a s ,« ® 1 + a\ R a Si g 1 + 10 a SiK a* - 
+ 10 - a S;S + a*- 10 a SiK + a SjS 10 a* ^ 



where {-, ■} denotes the anti-commutator. Comparing this with the generator of time translations 



H = cV = ^2 hujfi 



© (H a «,«' a L«J 0101 + 10 H a «,«' a i,J 01 + 1010 H a s,«> a L}) 



(58) 

(59) 
(60) 

(61) 
(62) 

(63) 
(64) 

(65) 



(66) 



(67) 



we can see that there is a relation between Tt and H but the latter contains terms describing interactions between the 
oscillators. The contribution from these interactions vanishes on vacuum states. Below, when we introduce the notion 
of a generalized coherent state, we will be able to relate averages of H and H. In a similar way one can introduce the 
"Pointing operator" 
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P = y d 3 x£{t,x)x B(t,x) = ^Y^ M ( a iz a sM + a s Mal^. (68) 

Its relation to the generator of 3-translations V is similar to this between TL and H. 

In the above construction the only element which is beyond a simple transition to many oscillators is the choice of 
c„. For different choices of these constants we obtain different representations of non-CCR and therefore also different 
quantization schemes. Several different ways of reasoning lead to c„ = 1/y/n as we shall also see in the next sections. 



V. SOME PARTICULAR STATES 



We assume that all the multi-oscillator states are symmetric with respect to permutations of the oscillators. 



A. Generalized coherent states 



For general c„ an eigenstate of a s ^ corresponding to the eigenvalue a s _z is of the form 

|a s ,«) = fi(s, k)\s, it, a S)i j/ci) 

©/2(s, k)\s, it, a s ^/ (2c 2 )) \s, it, a s ^/ (2c 2 )) 

®h{s, it)\s, it, a Si g/(3c 3 ))|s, it, a s ,g/ (3c 3 ))|s, it, a s ,g/ (3c 3 )) 



(69) 



where 

|s,K,a SiS ) = |s,K>|a SiK >, (70) 

S^Li l/n( s 5^)| 2 = 1, an d o|a S) g) = a s ,g|a s .g). What is interesting not all /„ have to be nonvanishing. 
The average "energies" of the field in the above eigenstate are 



x 1 I 

{cy. s .pi\n\a s .Fi) = huj R \a s ^\ 2 ^ — j\fn{s,K)\ 2 + -Tilor ^ n\f n (s, it)\ 2 (71) 

n=l nC " n=l 

and 

_^ 00 

(a a>it \H.\a a>it ) = hwjt\a s>it \ 2 + -fu^z ^ nc 2 n \f n {s, k)\ 2 . (72) 



2 „ , 



The two averages will differ only by the value of the vacuum contribution if c n — 1/ y/n which leads us back to the 
above mentioned choice of c„. With this choice and taking the general combination of coherent states 

i*> = X! $s .«i as .«) ( 73 ) 



we find 



1 00 



2 

s.K n=l 



and 



(*|H|*> = E &wa|*.,a| 2 Ka| 2 + J E ^a|<M 2 - ( 75 ) 



One may wonder, then, what is the more natural choice of the free-field Hamiltonian: H describing interacting 
oscillators, or H describing the noninteracting ones? The coherent-state average of H does not depend on the average 
number of the oscillators and naturally includes the process of energy exchange between different oscillators. Moreover, 
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both H and P are denned in the standard way in terms of the multi-oscillator non-CCR algebra. With this choice of 
free dynamics we find 

e lHt/h a s ^e- mt/h = e'^^a^ (76) 

as opposed to the standard formula 

e mt / h a s ^e- mt / h = e'^a^. (77) 

The latter choice is simpler because it leads to the standard form of the interaction-picture Hamiltonian and therefore 
will be the basis of our non-canonical quantum optics. The version based on H is a subject of ongoing study. 

B. Vacuum 

Similarly to the one-oscillator case the traditional notion of a vacuum state is replaced in our formalism by a vacuum 
subs-pace consisting of all the vectors annihilated by all a s ^. Their general form is 



|o> = £°S,ol*.*.°> 

S,K 

®E O S S2 ^,K 2 A0l S l^l'°)l S 2,«2,0) 



£ ,L2, S 3^ 1 ,KA 2 ,KA3,0,0,OI S 1:«1:°)I S 2,K2,0)|S3,K3,0) 

© . . . (78) 

It seems that there is no reason for introducing the standard unique "vacuum state" understood as the cyclic vector 
of the GNS construction. 

As an example of a vacuum state consider 

|0) = VpTIO) 

®^pl\0)\0) 
®s/pi\0)\0)\0) 

© . . . (79) 
The average energy of the free-field vacuum state is therefore 

oo 

H = (0|W|0) = ^2n Pn (0\H\0) =n~H (80) 

71=1 

where n and H are, respectively, the average number of oscillators and the average energy of a single oscillator. For 
the sake of completeness let us note that 

oo oo 

H= (0|H|0) = J2 c l n Pn(0\H\0) =Hj2<Z n Pn ( 81 ) 

n— 1 n— 1 

For c n — 1 we obtain H = H; for c n = l/\/n H = H the latter being independent of the number of oscillators. In 
both cases no problem with infinite vacuum energy is found. Obviously, one can contemplate also other vacua, say, 
in entangled or mixed states. 

C. Multi-photon states 

Assume |0) is a vacuum state. The non-canonical algebra (non-CCR) differes from the canonical one in the 
commutator 
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[ax,a{] = l x (82) 

where for any A, A' 

[l v ,a A ] = [l v ,4]=0 (83) 

where A stands for (s,k). A normalized state describing a collection of photons is defined in analogy to the standard 
formalism by 



1 



n 1 \n 2 \...n N \(0\lZ---ntlZ\° 



(a{ N r» . . . (a{ 2 r(alr\0) =: \n Xl ,n Xa , . . . ,n Ajv ) (84) 



States corresponding to the same A but different n's, or to the same n but different A's, are orthogonal. As a 
consequence a non-canonical vacuum average of any product of non-canonical creation and annihilation operators 
vanishes if and only if an analogous expression formulated in terms of the canonical objects does. This property is a 
consequence of three facts which hold true in both formalisms: (a) annihilation operators annihilate vacuum states, 
(b) creation operators are obtained by a Hermitian conjugate of the annihilation operators, and (c) the RHS of a 
commutator of creation and annihilation operators commutes with all creation and annihilation operators. 



VI. PERTURBATION THEORY 

It is essential that, similarly to the one-oscillator formalism, the free Hamiltonian (defined simply as a generator of 
time translations) generates the standard dynamics of annihilation operators: 

e mt/h a s ^ e - mt/h = e-^Oa.a. (85) 

Accordingly, the form of the interaction-picture Hamiltonian will by the same as in the standard theory. This would 
not be quite the same if we have chosen H in the role of the free Hamiltonian (an option which, nevertheless, should 
be investigated). In what follows we start with H = H n + V, where 

H = H A +H (86) 

V = --A(x)-p 
m 



m V 2ujfiV 

S,K ' 

In the interaction picture we get 

V(t) = --A(t,x)-p(t) 
m 



~ E J^V ( a ^e- i{ ^- R - S) e s ,, ■ Pit) + al^-^el, ■ p(tj) , (88) 



h_ 

S,K I 

and 

p(t) = e lHAt pe' lHAt . (89) 

Since we are purposefully neglecting the "A 2 " term in the Hamiltonian, one should restrict the analysis to the dipole 
approximation and therefore it is justified to set x = 0: 

= £(<**&(*) + «]$(*)). (90) 

A 

The operators 



are identical to those from the standard formalism and act only on atomic degrees of freedom (i.e. commute with 
a„s). 
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A. Spontaneous decay of an excited state 



The first problem we shall treat in the non-canonical way is a lifetime of an excited atomic state. The problem, as 
we shall see, is of particular importance for the physical interpretation of the non-canonical formalism. 

Assume that at t = the atom-field system is described by the state |\&(0)) = |0, A). The amplitude that the atom 
remains in the excited state is 

<(M|*(t)) = i 

+ dtzjjdh £<0|a A2 ajjO><A|& 2 (t 2 )£ Al (ti)| A) 

A1A2 

+ 7^1 / dt 4 [ 4 dt 3 [ 3 dt 2 [ 2 dh ]T (0|a A4 4 3 a A2 4jO)(A|ff A4 (t 4 )^ (ts)^^)^,^)!^) 
(in) J J J Q Jo , . , , 



Ai A 2 A3 A4 



+ ?4y4 / dU f* dt 3 f 3 dt 2 f 2 dh ( l a A 4 aA 3 ai 2 <|0)(A|.9A 4 (t4).9A3(i3).9i 2 (t 2 ).9i 1 (ti)|A) 

\ in ) Jo JO JO JO xxx\ " 

Ai A 2 A3 A4 

+ ... 
1 

+ / dt 2 [ 2 dh ^(OlaA.aijO)^!^^).^^)^)^^ 

^ ^ J° A 4 A 2 

+ ?4y4 / dt 4 [ 4 dt 3 [ 3 dt 2 [ 2 dh Yl (Ol^aLaA.aijO)^!^^)^^)^^)^^)^) 
l^J Jo Jo Jo JO \ \ a \ 



Ai A 2 A3 A4 



v T^OlOl 

x A A 4 A 3 A 2 Ai 



/ dt 4 [* dt 3 [ 3 dt 2 [ 2 dh V (6|a A4 aA 3 a A2 a Ai |6)(A|5 A4 (i 4 )5A 3 (t3)5i 2 (t2)5i 1 (ii)|^) 
Jo Jo JO JO \ \ \ \ 



(^) 4 Jo Jo 

yOOll 

X A A 4 A 3 A 2 Ai 



Ai A2 A3 A. 



+ . . . (92) 

In the above perturbative expansion we have explicitly shown all the nonvanishing terms up to the fifth order of 
perturbation theory. Here (6|a A4 a A3 4 2 4i |6) etc. are the expressions one would have obtained in ordinary canonical 
formalism and 

a a 4 a 3 a 2 a, - ~m — : — ~r~rrir 



L A 4 A 3 A 2 Ai „ „i ^ + I*, 

(0|aA 4 aA3aA 2 a Al |0) 



etc. Such expressions arc well defined since whenever their denominators vanish the whole term it containing van- 
ishes as well. This fact is of crucial importance and shows that the perturbative expansions in both canonical and 
noncanonical frameworks contain terms of exactly the same type but differing by the numerical factors X ". 

Let us note that in the above calculation we have not used the explicit realization of the non-canonical algebra but 
only the algebra itself. At such a general level both the canonincal and non-canonical theories can be regarded as 
particular cases of a more general theory characterized by the algebra 

[ax,a{,}=S X x'Ix, (94) 

[ax,a x >}=0, (95) 

[4,a A ,] =0, (96) 

[a{,lx'}=0, (97) 

[ax,h>}=0. (98) 

The canonical choice, based on oscillators with classical parameter A, is J A = 1; the choice based on oscillators with 
eigenvalue A is J A = 1 A . 

To proceed further and get more information as to the physical meaning of the non-canonical dynamics we have to 
make the analysis less general. First of all let us stick to the particular choice of 1 A in terms of POV measures we 
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have introduced earlier and assume that J2\ 1a = 1. Second, let us take the vacuum state in the form (|7^). Under 



such assumptions we can explicitly compute the factors X' 



TL-0101 

^ AAA' A' 



V0101 
^AAAA 



,-' . t = (0|1 A |0) - \O x \ 



,)\a x a{\6) 

= (Ol^lvlO) = $>(l - ~)W\0^ = (l - (l/n))M\0„\> 
(0|a A a A a v a A ,|0) ^ V nJ K '> 

(0\a\a\ a\a\ |0) , , , . 



,U.*t|o) 



<l/n))|0 A | 4 + (l/n)|0 A | 



AAAA 



vOOll 
^AA'AA' 



)|a A a A a A a A 

^sHW-(0|iJ|0)-(i-(i/»>)|O» 

(0|a A a A a A a A |0) v 7 

\L - \ t .t .ft/ = <0|1a1a-|0) = (l - <l/»>) |Oa 
(0|a A a A /a A a A ,|0) v 7 



<l/n)|0> 



2 |Ov| 



(99) 
(100) 
(101) 
(102) 
(103) 



What is interesting (and very characteristic) all these factors are smaller than 1 (this follows trivially from J2\ \0\\ 2 — 
1). An analysis of higher order terms shows that this is a generic property of the non-canonical perturbation theory. 
The n occuring in the average (1/n) is the number-of-oscillators operator. For realistic vacua one may expect the 
average number of oscillators to be large and therefore (1/n) w 0. Taking a more general vacuum state we arrive at 
(l/n A ) instead of (1/n), which means that the influence of the vacuum may vary from frequency to frequency (i.e. 
from point to point in space). 

Now, if we ignore the corrections coming from (1/n) we can see that the non-canonical perturbative expansion of the 
amplitude is the same we would have obtained by using the standard theory but with g\ regularized by g\ — ► 0\g\. 
As a consequence there exists a natural cut-off in the theory which follows only from the fact that the vacuum wave 
function is square-integrable and therefore O s ^ — > for \k\ — > oo. It is quite remarkable that the same mechanism 
that eliminated the infinite vacuum energy plays a similar role in the other parts of the theory. As we shall see 
shortly the actual role of the vacuum can be analyzed only in a fully relativistic setting since then the charge and 
mass renormalization come into play. 

However, for the sake of concreteness and to make some rough estimates of the effects involved let us take the 
trivial example where the vacuum amplitudes are constant, say, O s ji = C for all < to max and zero otherwise. 
The dynamics of the amplitude is then (up to (1/n) ss 0) equivalent to the standard one with the cut-off at w max 
and the coupling constant e/m replaced by Ce/m. This implies that e and m have to be treated as bare parameters 
and the experimantal value is e ex /m ex = Ce/m. With this observation in mind we can discuss non-canonically other 
perturbative effects which are widely believed to be a consequence of the standard canonical quantization. Below 
we find it useful to make the bare parameter e/m explicit in calculations and for this reason will use the notation 
gx = (e/m) /a. 



B. Spontaneous emission of N identical photons in iV-th order perturbation theory 



Assume the atomic transition is \A) — > \B). Up to the N-th order perturbative correction we get 

/ft P^N ft, 

o 

(^y%l\f4) X (-levantcanonicalformula) 

N 



(N x ,B\*(t)) = (^-^ fdt N ... r dt 1 (N x \a{ N \0)(B\fi(t N )...fi(t 1 )\A) 



(-) 



[0|lj^|0) x (relevant canonical formula) 



(104) 



where the "hatted" expressions are those from the canonical theory. As we can see the task is reduced to computing 
(0 j ljY |0). The general formula, valid for any N, is somewhat complicated and not very illuminating. The cases N = 1 
and N = 2 we have already met. Making the simplifying choice of a very "flat" distribution of the vacuum modes 
(i.e. O s ^ = C below some threshold) we get 
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(l x ,B\*(t)) = C{t x ,B\9(t)) (105) 
(2 A , B|tt(t)> = sjc* + <l/n>(C* - C^)(2 Al B|*(t)>. (106) 

Identifying Ce/m with e cx /m cx we obtain the standard quantum-optics results provided 

y/c* + (l/n)(C 2 -C*) « C 2 . (107) 

For N = 3 

oo 1 

(0|1 A |0) = ^ — (nC 2 + 3n(n - 1)C 4 + (n 3 - 3?i 2 + 2n)C* 6 )p„ 

71 = 1 71 

= C 6 + 3C 4 (l - C 2 )(l/n) + C 2 (l - 3C 2 + 2C 4 )(l/n 2 ). (108) 
As before the result becomes the standard one if the approximation 

^C 6 + 3C 4 (l - C 2 )(l/n) + C 2 (l - 3C 2 + 2C 4 )(l/n 2 ) « C* 3 (109) 

is justified. Let us note that in the standard cananonical quantum optics one considers vacuum consisting of an 
infinite number of oscillators and such subtleties are trivially ignored. 

C. Spontaneous emission of two different photons 

By the same argument as before the proportionality factor we need to estimate in second-order perturbation theory 

is 

(l x ,l x ,\a{a{,\0) 

<i A ,i A 'l44'|6) 



= v/(0|l A lv|0) - J I - (l/n)C 2 « C 2 (110) 



A A' 

so the result agrees with the canonical one. It is also in perfect agreement with the explicit calculations given in |J 

D. Stimulated emission 

The last example we will discuss is the first order calculation of the transition amplitude |2V A , A) — > |(JV + 1) A , B) 
The appropriate proportionality coeeficient is 



provided the values of (l/n), (l/n 2 ), etc. are sufficiently small. This is the correct result since C gets absorbed into 
the renormalized coupling constant. 

VII. BLACKBODY RADIATION 

The final test of the new formalism we want to perform is the problem of blackbody radiation. Planck's famous 
formula Q 

j . 7i i*u 3 h , . 

= ^2 C 3 e 0Hu _ i = N »> ( 112 ) 

where is the average number of excitations of an oscillator in inverse temperature /3, is one of the first great 
sucesses of quantum radiation theory and marks the beginning of quantum mechanics. Contemporary measurements 
of q{ui) |l^,|l3| performed by means of COBE (Cosmic Background Explorer) are in a very good agreement with 
the Planck law. The data have been carefully analyzed in the context of nonextensive statistics HJlql in search of 
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possible deviations from extensivity. The result that comes out systematically is \q — 1| < 10~ 4 where q is the Tsallis 
parameter. The case q = 1 corresponds to the exact Planck formula. If there are any corrections whatever, they must 
be quite small. 

The derivation of the formula given by Einstein is based on the properties of spontaneous and stimulated emissions. 
As we have seen above there may occur differences with respect to the standard formalism but under reasonable 
assumptions they may be expected to be small. 

Below we follow another standard route which consists basically of two steps. First, one counts the number of 
different wave vectors k such that c\k\ € [u),oj + Aw]. Second, one associates with each such a vector an oscillator and 
counts the average number of its excitations assuming the Boltzmann-Gibbs probability distribution at temperature 
T and chemical potential /i = 0. The latter assumption is justified by the fact that the number of excitations of the 
electromagnetic field is not conserved in atom-light interactions. 

In the new model the situation is slightly different since there exists an additional conserved quantum number: 
The number of oscillators . As we have seen in previous calculations the Hamiltonian is block-diagonal with respect 
to © but changes the number of excitations in each ./V-oscillator subspace of the direct sum. The state vectors at 
the multi-oscillator level are symmetric with respect to permutations of the oscillators and therefore the oscillators 
themselves have to be regarded as bosons whose number is conserved and their chemical potential is /i ^ 0. However, 
their excitations should be regarded as bosons with vanishing chemical potential. 

The eigenvalues of Ti 

E m ,n = mhiu(n + l^J . (113) 

corresponding to the oscillator whose frequency is u> are parametrized by two natural numbers: m (the number of 
oscillators) and n (the number of excitations). Assuming the standard Boltzmann-Gibbs statistics we obtain the 
probabilities 

Pm,n = Z- 1 e- fJ[mhuj( - n+ ? ) - rn > l] (114) 

where 

3 — flmhuj 



e /?m( M +W2) e (115) 



rn—1 



The Lambert series |16 



^ X* 

1^ °™— 



(116) 



m— 1 



is convergent for any x if X)m=j a ™ ^ s convergent. Otherwise (116) converges for exactly those x for which the power 



series Y,Z=i a mX m does - In (HH) »™ = e^" 1 ^^^ and £~ =1 a m < oo if \i + Tilu/2 < 0. If [i + Tiuj/2 > we still 
have convergence of (115) as long as Ylm=i e - ' 3 " 1 ^''" - ^] < oo. The upper limit imposed on \x by the finiteness of Z 
is therefore /i < ^huj. In what follows we assume that fi is w-independent and therefore fi < 0. 
The appropriate average number of excitations is 



oo oo 



z -i J2 J2 mne~ p[mhuj(n+ i ) - mfl] (117) 



m— 1 n— 

and the Planck formula is replaced by 



0newM = —^UJ^n^. (118) 

It is easy to show that £> n ew(w) tends to the Planck distribution with /i — > — oo. To see this consider a more general 
series 

oo oo 

Z- 1 J2 mnq m e~ f}mhuj( - n+ i ) (119) 
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where Z is the normalization factor and X)m=i 9m < oo. If Qi — 1 an d 9m = for m > 1 then (119) is just the exact 



Planckian formula. Factoring out e 
q m = e-ZSHC™- 1 ) for m > 1. For \li\ 



-P\iA i n both the numerator and the denominator of n w we obtain q\ = 1 and 
oo all <7 m , for m > 1, vanish and the limiting distribution is Planckian. 
This proves that an experimental agreement with the ordinary Planck's g(co) cannot rule out our modification 
but can, at most, set a lower bound on an admissible value of \fi\. However, assuming that li has some finite 
and fixed value it should be in principle measurable. The plots show that the modifications become visible around 
fi w — 3ksT. Assuming that the chemical potential is temperature independent, say fi — — fcsTo, we obtain a kind 
of critical temperature T cr iti ca i ~ To/3 above which the ratio /i/(fcsT) is small enough to make the modifications of 
the distribution observable. For T < Tcriticai the distribution should be given by the Planck law; for T > T cr iti ca i the 
distribution should approach the ji = distribution, i.e. this would be a Planck-type curve but with the maximum 
lowered and shifted towards higher energies. 

Fig. 1 shows the plots of g nev/ (u>) for fi — (lower dotted), ll — — 0.8fcsT (upper dotted), and fi = — lOfc^T (solid). 
The thick dashed curve is the Planck distribution. The curve obtained for ll — — lOfc^T is indistinguishable from the 
Planck distribution. The plot does not change if one takes ll < — lOfc^T and differences are not visible even if one 
plots the distributions in logarithmic scales (not shown here) . This is a numerical proof that the distribution we have 
obtained on the basis of the modified quantization tends very quickly to the Planck one as fi — > — oo. It is instructive 
to compare the modification we have predicted with those arising from nonextensive statistics. The two thin dashed 
lines represent Tsallis distributions resulting from nonextensive formalism for q = 0.95 (lower) and q = 1.05 (upper). 
The modifications we have derived are therefore qualitatively different from those resulting from Tsallis statistics. 




2 4 6 8 10 

FIG. 1. Q ncv (uj) for jj, = (lower dotted), \i = — 0.8/cgT (upper dotted), and /i = — 10fc_eT (solid). The energy range 
is O.OlfcsT < hcj < lOfcsT. The thick dashed curve is the Planck distribution. The curve obtained for fj, = — lOfcsT is 
indistinguishable from the Planck distribution. The two thin dashed lines represent Tsallis distributions resulting from the 
Tsallis formalism for q = 0.95 (lower) and q = 1.05 (upper). Since g ncv < g at least in the neighborhood of the maximum, 
the new distribution has to be compared with q < 1 statistics. The curves are qualitatively different. In particular, all q < 1 
distributions require an energy cut-ofT which does not occur for g nC w 
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FIG. 2. encw(^) for -10 < -T /T < 0. The cut through T /T = 10 is practically indistinguishable from the Planck 
distribution. 
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FIG. 3. Contour plot of Q uevl (u>) for 
clearly seen. 




6 8 10 

-10 < —Tq/T < 0. The fast convergence to Planckian q(oj) (as Tq/T increases) is 



VIII. CONCLUSIONS 



"A theory that is as spectacularly successful as quantum electrodynamics has to be more or less correct, although 
we may not be formulating it in just the right way" ]l0| . The above quotation from Weinberg could serve as a motto 
opening our paper. The main idea we have tried to advocate was that the standard canonical quantization procedure 
is, in certain sense, too classical to be good. 

The reasons for such a choice of quantization could be both historical and sociological and may be rooted in the fact 
that the idea of quantizing the field was formulated before the real development of modern quantum mechanics. In 
oscillations of a simple pendulum it may be justified to treat u> as an external parameter defining the system (via, say, 
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the length of the pendulum). But oscillations of the electromagnetic field do not seem to have such a "mechanical" 
origin and it is more natural to think of the spectrum of frequencies as eigenvalues of some Hamiltonian. That is 
exactly what happens with other quantum wave equations. 

We have defined the quantum electromagnetic field as an oscillator that can exist in a superposition of different 
frequencies (or, rather, wave vectors). This should not be confused with the classical superpositions of frequencies 
created by, say, a guitar string. The superpositions we have in mind dissappear at the classical level. 

Once one accepts this viewpoint it becomes clear how to quantize the field at the level of a single oscillator. We 
do not need many oscillators to perform the field quantization. But there is no reason to believe that all the possible 
fields can be described by the same single oscillator. And even more: We know that the structure of the one-oscillator 
Hilbert space is not rich enough to describe multi-particle entangled states and there is no doubt that such states 
are physical. The next step, performed already after the quantization, is to consider fields consisting of 1, 2, 3 and 
more oscillators, and even existing in superpositions of different numbers of them. The resulting structure is in many 
respects analogous to the Fock space so that the procedure can be (although somewhat misleadingly) referred to as 
"second quantization" . What is essential we do not need the vacuum state understood as the unique cyclic vector of 
the GNS construction. 

On the other hand, there exist vacuum states. These are all the states describing ground states of the oscillators. 
They correspond to concrete finite average values of energy. A general vacuum state is therefore a superposition of 
different eigenstates of a free Hamiltonian and is not, in itself, an eigenstate of the Hamiltonian. 

One technical assumption we have made is the resolution of identity property of the non-CCR algebra. This 
assumption is clearly satisfied at the one-oscillator level. The remaining assumptions are standard. The system is 
described by laws of ordinary quantum mechanics so that to compute concrete problems we can use standard methods. 
Perturbation theory leads to structures we know from the standard Feynmann diagrams. The blackbody radiation is 
calculated by means of the standard Boltzmann-Gibbs statistics. 

There are still several unexplored possibilities. To give an example, we have seen in explicit calculations that differ- 
ences between the canonical and non-canonical formalisms consist of the factors occuring in perturbative expansions 
of amplitudes. These factors explicitly depend on the choice of the RHS of the non-CCR algebra. As such, they 
point into possible experiments testing directly the algebra of canonical commutation relations. The meaning of such 
tests is, at the present stage, obscured by the lack of proper understanding of the role of renormalisation. Actually, 
one should not expect here precise results since even the nonrelativistic canonical quantum electrodynamics is non- 
renormalizable. The simple illustration we have used, namely the one with the cut-off and flat vacuum, shows that 
the new theory is not that far from the canonical one as one might expect. 

Let us close these remarks with another quotation: "Present quantum electrodynamics contains many very impor- 
tant 'elements of truth', but also some clear 'elements of nonsense'. Because of the divergences and ambiguities, there 
is general agreement that a rather deep modification of the theory is needed, but in some forty years of theoretical 
work, nobody has seen how to disentengle the truth from the nonsense. In such a situation, one needs more exper- 
imental evidence, but during that same forty years we have found no clues from the laboratory as to what specific 
features of QED might be modified. Even worse, in the absence of any alternative theory whose predictions differ 
from those of QED in known ways, we have no criterion telling us which experiments would be relevant ones to try. 

It seems useful, then, to examine the various disturbing features of QED, which give rise to mathematical or concep- 
tual difficulties, to ask whether present empirical evidence demands their presence, and to explore the consequences 
of the modified (although perhaps rather crude and incomplete) theories in which these features are removed. Any 
difference between the predictions of QED and some alternative theory, corresponds to an experiment which might 
distinguish between them; if it appears untried but feasible, then we have the opportunity to subject QED to a new 
test in which we know just what to look for, and which we would be very unlikely to think of without the alternative 
theory. For this purpose, the alternative theory need not be worked out as completely as QED; it is sufficient if we 
know in what way their predictions will differ in the area of interest. Nor does the alternative theory need to be free 
of defects in all other respects; for if experiment should show that it contains just a single 'element of truth' that is 
not in QED, then the alternative theory will have served its purpose; we would have the long-missing clue showing in 
what way QED must be modified, and electrodynamics (and, I suspect, much more of theoretical physics along with 
it) could get moving again " (I). 
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